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SUMMARY 


A semiempirical relationship is presented which describes the extent of 
interaction between constituents in single-phase binary alloy systems having 
planar, cylindrical, or spherical interfaces. This relationship makes possi- 
ble a quick estimate of the extent of interaction without lengthy numerical 
calculations. It includes two parameters which are functions of mean concen- 
tration and interface geometry. Experimental data for the copper -nickel sys- 
tem are included to demonstrate the usefulness of this relationship. 


INTRODUCTION 

In the analysis of diffusion controlled processes (e.g. , the diffusion of 
planar protective coatings into substrate materials, the filament-matrix inter- 
action in metal matrix composites, or the homogenization of powder compacts) , 
the ability to predict the degree of interaction between different components, 
given some exposure conditions, is needed. In the past, the approach has been 
to solve the diffusion equation with the appropriate initial and boundary con- 
ditions (refs. 1 and 2). The solution to the diffusion equation is usually 
obtained by numerical methods such as finite-difference or infinite-series 
solutions except in limited cases where analytical solutions exist. However, 
the numerical methods involve extensive calculations, and a separate set of 
calculations is required for each composition and interfacial geometry of the 
material system being considered. 

The purpose of this paper is to provide a single semiempirical relation- 
ship which describes the extent of interaction between the constituents in 
single-phase binary alloy systems having planar, cylindrical, or spherical 
interfaces. This relationship makes possible a quick estimate of the extent 
of interaction without lengthy numerical calculations. It includes two 
parameters which are functions of mean concentration and interface geometry. 
Experimental data for the copper-nickel system are included to demonstrate 
the usefulness of this relationship. 


ANALYSIS 

Figure 1 shows the nomenclature and schematic of the geometric models for 
planar, cylindrical, and spherical interfaces. The average concentration H 
of element B in the original B-rich region is a measure of the degree of homog- 
enization at any time t. For the planar, cylindrical, and spherical geome- 
tries (denoted by the subscripts p, c, and s, respectively) , H is given by 
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where Cg is the volume fraction of element B at distance r, Rg is the 
thickness or radius of the B-rich region, and R T is the total thickness or 
radius of the sample. (See fig. 1.) The degree of homogenization can also be 
expressed as the ratio of the amount of material transferred across the bound- 
ary at time t to the amount of material at t = 00 (ref. 2) . The relation- 
ship between this ratio F and H is given by 


H = 1 - F(1 - C B ) (4) 

where Cg is the average volume fraction of element B in the specimen. The 
value of F varies from zero to unity whereas H varies from unity to Cg. 

If the diffusion coefficient D is assumed to be independent of concen- 
tration, the equation governing homogenization in single-phase binary alloy 
systems is given by 


9C B = _D_ l_( m 3Cb ) 
9t jjn 3r V dr J 
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where m is 0, 1, and 2 for planar, cylindrical, and spherical geometries, 
respectively. The initial and boundary conditions to be satisfied are 
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Solutions to equations (5) and (6) were determined for each geometry for 
different values of Cg by using a finite-difference technique (ref. 3) . 

Values of H were plotted as a function of x for each geometry and Cg, 
where T is a dimensionless parameter and is equal to the normalized diffusion 
time Dt/Rg 2 . Several different equations were fitted empirically to the 
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plotted values to determine the simplest expression that best described 
all data. 

For the special case of planar geometry with c B + 0 , a simple closed- 
form solution exists for equations (5) and (6). In this case, the concentra- 
tion of component B (ref. 4) is given by 
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where erf and erfc are the error function and the error function compli- 
ment, respectively. The equation in terms of erfc is in a convenient form 
for the subsequent integration required in developing an expression for Hp, 
which is the homogenization parameter given by equation (1) . By making use of 
the relation 


1 -T - 2 

erfc ( T) ) df| = — e - r erfc (r) 


equation (1) reduces to 



For this case, the empirical results were compared with those from the analyti- 
cal solutions (eq. (9)) as well as with those from the finite-difference 
solutions. 


RESULTS AND DISCUSSION 

The following empirical equation gave results which were in good agreement 
with those obtained from the finite-difference analysis: 



H = C B + (1 - C B ) exp(-pt q ) 


( 10 ) 
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Parameters p and q were evaluated for each geometry and value of C B that 
were considered. 

Typical results for planar geometry with C B = 0.25 are shown in fig- 
ures 2 and 3, where the curves of H as a function of T which were determined 
from equation (10) are compared with those from the finite-difference analysis. 
Curves are shown for different values of q (fig. 2) and p (fig. 3) to deter- 
mine the sensitivity of H to these parameters. As can be seen, p = 1.15 
and q = 0.65 gave results nearly identical to those from the finite-difference 
analysis. Similar calculations were carried out for cylindrical and spherical 
geometries, and in all cases equation (10) satisfactorily described the finite- 
difference results. 

Table I shows the p and q values for all geometries considered. The 
quantity Ah is the absolute difference between a finite-difference value 
and the corresponding value computed with equation (10) . As can be seen from 
table I, the values for AH max and AH av (where the subscripts max and av 
denote maximum and average, respectively) are less than or equal to 0.033 and 
0.012, respectively, with the exception of planar geometry with C B = 0.01, 
for which the values are as high as 0.070 and 0.025, respectively. The ana- 
lytical solution (eq. (9)) could be used if greater accuracy were desired. 

Figure 4 shows a plot of p and q as a function of C B for each of the 
three geometries. Parameter p increases with increasing C B . Also, for a 
given concentration, the values of p for cylindrical and spherical geometries 
are, respectively, about 2 and 3 times as great as that for the planar geometry. 
However, parameter q has a maximum for each geometry. The values of C B at 
the maxima are approximately 0.5, 0.3, and 0.22 for planar, cylindrical, and 
spherical geometries, respectively. For planar geometry, the curve of q as 
a function of C B is symmetrical about C B = 0.5. This value of C B corre- 
sponds to diffusion conditions where the thicknesses of the A and B regions 
are equal (Rj, = 2R B ) . The asymmetry in the cylindrical and spherical curves 
is due to the complexities of the interface geometry. 

Carpenter, Tenney, and Houska (ref. 5) have electroplated 8 ym of nickel 
onto one face of a (111) oriented copper single-crystal disk, annealed the 
specimen at 1173 K for 3/4, 2, and 5 hr, measured diffracted X-ray intensity 
after each anneal, and determined concentration profiles from the X-ray data 
by an absorption technique. Because the plating thickness (R B = 8 ym) is 
negligible when compared with the total thickness R T , the degree of homogeni- 
zation in this specimen can be calculated by using the values of p and q 
that correspond to C B = 0 (i.e., p = 0.55 and q = 0.50). The solid line 
in figure 5 denotes the variation of H with T which was generated by using 
these p and q values in equation (10) . For this planar geometry, the 
empirical curve is in reasonable agreement with the finite-difference and exact 
solutions as well as with the experimental data points. The diffusion coef- 
ficient D was taken as 0.2 x 10 “-^ m ^/s, which corresponds to the nickel-rich 
concentration range (ref. 6) . A slightly lower value of D would be more 
appropriate because 0.2 x lO - -^ m^/s corresponds to polycrystalline nickel 
whereas the present specimen consisted of nickel epitaxially deposited on 
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single-crystal copper. Renormalizing T with the lower D value would move 
the experimental points toward lower T values, making the agreement better. 

Heckel (ref. 2) prepared compacted mixtures of elemental copper and nickel 
powders with different average compositions, annealed them in the 1223 to 1323 K 
temperature range for various lengths of time, and analyzed them for the degree 
of homogenization F. His homogenization results F for a sample with 
C B = 0.11 were converted to H by using equation (4) . The results are shown 
in figure 6 and can be compared with the finite-difference and equation (10) 
results if a spherical geometry is assumed. The values of p (3.8) and 
q (0.66) used in equation (10) were obtained from figure 4 for C B = 0.11. 

As can be seen, the present empirical solution is in excellent agreement with 
the finite-difference solution and is in reasonable agreement with the experi- 
mental homogenization data on powder compacts. 

Figure 7 shows the predicted homogenization values for specimens with 
cylindrical geometry and C B = 0.01, 0.50, and 0.75. No experimental data for 
long diffusion time could be found in the literature for the cylindrical geom- 
etry case although it would be of importance in the analysis of interaction in 
fiber-reinforced composites. (Herring and Tenney (ref. 7) have investigated 
diffusion between Ni filaments and a Cu matrix. However, the extent of inter- 
action in their studies is not large enough to provide a meaningful check of 
the results shown in fig. 6.) The finite-difference solutions for C B = 0.50 
and 0.75 are identical to the present empirical solution and are not shown for 
T < 10 -1 to avoid confusion. As can be seen, the present empirical solution 
is in good agreement with the finite-difference solution for all the values 
of C B considered. 


CONCLUSIONS 

A semiempirical relationship has been developed which describes the extent 
of interaction between constituents in single-phase binary alloy systems with 
planar, cylindrical, or spherical interfaces. The relationship contains two 
parameters which are functions of mean concentration and interface geometry of 
the couple. The empirical solution is simple, easy to use, does not involve 
sequential calculations, and hence makes possible a quick estimate of the 
extent of interaction without lengthy calculations. Results obtained with 
the empirical relationship are in good agreement with results from a finite- 
difference analysis and with experimental results. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton; VA 23665 
October 4, 1978 
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TABLE I.- RESULTS OF CURVE FITTING FOR DIFFERENT 


C B 


AND INTERFACIAL GEOMETRY 


Geometry 

^B 

P 

q 

^ H max 

AH av 

Planar 

0.01 

0.61 

0.50 

0.070 

0.025 

Planar 

.10 

.80 

.55 

.028 

.012 

Planar 

.25 

1.15 

.65 

.011 

.004 

Planar 

.50 

2.25 

.69 

.014 

.007 

Planar 

.75 

4.50 

.63 

.005 

.001 

Planar 

.90 

8.98 

.56 

.004 

.002 

Planar 

.99 

30.30 

.43 

.001 

.000 

Cylindrical 

.01 

1.62 

.61 

.033 

.009 

Cylindrical 

.10 

2.14 

.65 

.012 

.004 

Cylindrical 

.25 

3.05 

.67 

.020 

.009 

Cyl indrical 

.50 

5.08 

.66 

.011 

.005 

Cylindrical 

.75 

8.42 

.58 

.005 

.001 

Cylindrical 

.90 

14.40 

.52 

.005 

.002 

Spherical 

.01 

2.95 

.64 

.013 

.003 

Spherical 

.10 

3.74 

.66 

.017 

.008 

Spherical 

.25 

5.13 

.67 

.020 

.009 

Spherical 

.50 

8. 14 

.65 

.009 

.003 

Spherical 

.75 

13.19 

.59 

.007 

.003 




SPHERICAL 



C B = (Rg/y 3 


Figure 1.- Geometric models showing zero flux boundary conditions 
(J = 0) and average concentration C B . 
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T 


Figure 2.- Sensitivity of H to changes in parameter q for planar 
geometry with C B = 0.25 and p = 1.15. 



T 


Figure 3.- Sensitivity of H to changes in parameter p for planar 
geometry with C R = 0.25 and q = 0.65. 





H 



T 


Figure 5.- Experimental and predicted homogenization values for 
copper-nickel planar couple (C N ^ - 0) . 



Figure 6.- Experimental homogenization data for copper-nickel powder 
compact and the predicted values if a spherical geometry is 
assumed (C M - = 0.11). 
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